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1 Introduction 



The top quark pair production close to the threshold is an important process at a future 
International Linear Collider (ILC). It can be used to determine top quark properties, 
like the mass rrit and the width Ft, but also the strong coupling a.^ with high precision. 
This is in particular true for nit where an uncertainty below 100 MeV can be obtained 
from a threshold scan of the cross section [1]. 

The feasibility of such high-precision measurements requires a theory prediction of the 
total cross section a{e~^e~ — tt) with high accuracy (preferably 5a /a < 3%). Current 
estimates based on (partial) next-to-next-to-leading logarithmic (NNLL) order [2,3] and 
(partial) next-to-next-to-next-to-leading order (NNNLO) [4, 5] QCD corrections lead to 
an uncertainty of the order of 10%. 

In order to reach a theory goal of 5a/a < 3% it is necessary to include in the prediction 
next to the one-loop electroweak corrections, which are known since quite some time [6] 
(see also [7]), also higher order effects. The evaluation of 0{aas) corrections has been 
started in Ref. [8], where the two-loop mixed electroweak and QCD corrections to the 
matching coefficient of the vector current has been computed due to a Higgs or Z boson 
exchange in addition to a gluon. The current paper continues this enterprise and provides 
a result of 0{aas) for the two-loop vertex diagrams mediated by a boson and gluon 
exchangelll This completes the vertex corrections of order aus — a building block for 
the top quark production cross section. Assuming the (numerically well justified) power 
counting a ~ one can see that these corrections are formally of NNNLO. 

In order to complete the matching corrections of order aag also the two-loop box diagrams 
contributing to e^e~ ti have to be considered. Actually, only the proper combination of 
the box, vertex and self-energy contributions (the latter can, e.g., be found in Refs. [9,10]) 
forms a gauge independent set. 

The remainder of the paper is organized as follows: In the next Section we introduce our 
notation and derive the cross section formula for e^e~ ttX near the tt threshold. We 
present a general formula which includes all radiative corrections of the Standard Model 
(SM). In Section [3] we discuss some technical details of the two-loop computation and 
in Section m we concentrate on the 0{aas) corrections to the jti vertex and present our 
results. Section [5] contains our conclusions. Additional useful material concerning the 
one-loop expressions can be found in the Appendix. 

2 Threshold cross section 

The production cross section for the process e'^e~ tiX near threshold consists of 
helicity amplitudes and the hadronic part. The former correspond to the high-energy 

^Of course, in addition to the gauge boson also the corresponding Goldstone boson is taken into 
account. 
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production amplitude of a top quark pair, the latter to the QCD bound-state dynamics 
of the produced tt pairs exchanging gluons to form a resonance. The cross section can be 
cast in the form (for left-landed positrons and right-handed electrons), 

R{ele^-^tiX) = ^lm[{hn,v)''Hv + {hn,AfHA\, (1) 

where s is the square of the center-of-mass energy and -R(e^e^ tiX) is the cross section 
normalized to (j{e^e~ f^'^fJ'^) = (4'7ra^)/(3s). The first subscript of h refers to helicity 
of the electron, and the second one to the vector {Jy = ip'^^ip^oi axial- vector coupling 
(J^ = of the gauge bosons to the top quark current o For e^e^ in the initial 

state a similar expression is obtained by replacing R by L in Eq. ([T]). In the SM the 
helicity amplitudes for e^e^ and e^e^ are proportional to uie/Mw ~ 10~^ and are thus 
negligible. 

In the cross section formula ([1]) "Im" refers to those cuts which correspond to the tiX final 
statej^l This means that we have to select special cuts which correspond to the final state 
we are interested in. This requires a dedicated study incorporating the experimental setup. 
For the one-loop electroweak correction this treatment was performed in [11]. In this paper 
we will not pursue this problem further (see also the discussion in the Conclusions). 

The hadronic part is described by non-relativistic QCD (NRQCD) [12]. For our purpose 
it is sufficient to re-write the vector current Jy and the axial- vector current in terms of 
two-component NRQCD spinor fields ip, Xi which correspond to non-relativistic top and 
anti-top quarks, respectively. This yields the following NRQCD currents 

With the help of the NRQCD equation of motion for top and anti-top quarks (^/'V* D x = 
mtidolip'^a^x)) the 1/m^-suppressed vector current can be re-expressed in terms of 
Thus our matching relation between SM and NRQCD currents are given by 



Jy = e^^-*-o (^c, - ^z^o J , Ja = e^^™'^«c, j^, (3) 

with c^, = (i„ = Cq = 1 at tree level. The hadronic part is defined by the current correlation 
function 

Hf = Uxe^'^^\n\i:j'p\x)jmm {F = V,A), (4) 

1. 'J 



^The notation is basically adapted from Ref. [6], however, we added a second subscript F = V^Ato 
incorporate the axial- vector coupling of the Ztt vertex (see Eq. ([5])). 

■^From the theoretical point of view one has a pure ti final state up to NNLO in QCD. Starting from 
NNNLO one has to include the real emission of a gluon. Once the electroweak sector is considered, final 
states like {hW^)t needs to be included, where (hW^) has an invariant mass p^^^^ in a range \Vhw~''^t \ ^ 
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where E = ^/s — 2mt and is the NRQCD vacuum state. 

The evaluation of Hp requires to integrate out the low-energy modes of QCD, the soft, 
potential and ultrasoft gluons contained in NRQCD [13,14]. For the top quark system this 
can be done perturbatively. In a first step one integrates out the soft and potential gluons 
which results in the effective field theory Potential NRQCD [15,16]. The corresponding 
Lagrangian is known to NNNLO [17] To integrate out nonrelativistic top and anti-top 
quark fields the Rayleigh-Schrodinger perturbation theory can be applied as was initiated 
in Ref. [19] and performed to NNNLO for Hy and Ha in Refs. [4] and [20], respectively. 
Integrating out the ultrasoft gluon was completed recently in Ref. [5]. For the details of 
these steps we refer the reader to the original papers and references cited therein (see also 
Refs. [21-25]). 

In this paper we restrict ourselves to hard loop corrections to the production cross section, 
namely the corrections being parameterized as hj^p. The tree- level expressioiilfl of helicity 
amplitude hi^p is given by 



htv = QeQt + with Py 



s-Mz 



where the /?/ is the coupling of a fermion (/ = e, t) to the Z boson, is the sin of the 
weak-mixing (c^ = 1 — s^), and electric and iso-spin charges for top quark and electron 
are given by 

Qe = -1, Qt = 2/3, {n)t, = l/2, (r3)e, = -1/2, (T3)/^ = 0. (6) 

In the following the abbreviation = (T3) will be used. hj^A can be obtained by 
substituting /3y by (3^ = (/^^ — l3\}/2 in formula 

Let us now explain how hard loop corrections within SM can be incorporated into the 
helicity amplitude hi^p. To this end we organize the corrections as 

hi,A = hfl + hfl\ (7) 

where the h^{f and h!fip ■* (due to the jy^"^ "*) are the tree-level contributions, and h^j'p 
incorporate the contributions from radiative corrections (the superscript (z, j) denotes the 
electroweak- and QCD-loop order, respectively). As one can see from the expression of 

the axial- vector current ja-, Ha is suppressed by D jvn^ ~ E/mt. Thus one-loop QCD 



''The only missing constant in Ref. [17] is related to the three-loop static potential where recently the 
fermion corrections became available [18]. 

^We include the effect due to ^™ with cl^^°° = 1 into h^jij"' for convenience, see Eq. ([7]). 
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Figure 1: Two- line Mis. The solid and dashed lines correspond to propagators with mass 
rrit and M^, respectively. 



corrections to hj ^ correspond to NNNLO effects. Hard QCD corrections to the 'jtt and 
Zti vertices modify the matching coefficients Cy^a at loop level. We absorb these effects 
into helicity amplitudes and obtain (using hf(P "* = —h^/y E/{6mt)) 

/,(0.1) _ _ ^^^\ > tree > (0,1) _ (1) , tree /oN 

with Cv\ dv'^ being the z-loop contribution to the matching coefficients. For the purpose of 
this paper only the one-loop contribution clp is needed (see below for explicit expressions) . 
Let us mention that the two-loop QCD corrections have been evaluated in Refs. [26,27] 
and the three-loop corrections induced by a light quark loop in Ref. [28]. 



3 Technical details of the two-loop calculation 

Let us in this Section provide some technical details about the evaluation of the two- 
loop diagrams. They are generated with QGRAF [29] and further processed with q2e and 
exp [30,31]. The reduction of the integrals is performed with the program crusher [32] 
which implements the Laporta algorithm [33,34]. We arrive at 29 master integrals (MI) 
which are depicted in Figs. [IHH All diagrams occur with the propagators raised to power 
one. Note that there are two more Mis of type 3.10: one with a squared top quark 
propagator and one with a squared massless propagator. Similarly, an additional MI 
arises from type 3.11 with a squared massless propagator. 

We refrain from presenting the explicit results for all Mis in this paper but provide them 
in form of a Mathematica fil^ MIttewW.m using the conventions as defined in Eq. (136!) 
which corresponds to the one-loop tadpole integral. In all results presented in this file a 
factor (/i^/m^)^^ with c? = 4 — 2e has to be multiplied. 

Some Mis factorize into one-loop integrals or contain only one dimensionful scale. Most 
of these integrals are available in the literature and can, e.g., be found in Refs. [8,35-39]. 

As we will see in Section H] a rapid convergence is observed if one considers an expansion 
of the matching coefficient in the quantity z = M^/mf. For this reason we evaluate the 
two-scale Mis in this limit. A promising method is based on differential equations (see 

^The file is available from |http : //www-ttp . particle . uni-karlsruhe . de/Progdata/ttp08/ ttp08-43 . 
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Figure 2: Three-line Mis. The sohd, dashed and dotted hnes correspond to propagators 
with mass m^, and 0, respectively. Single external lines are on the mass shell with 
mass rrit whereas double external line have mass 2mt. The topology denoted by 3.10 
contains three master integrals, topology 3.11 contains two Mis. 




Ref. [40] for a recent review) which provide the expansion in an automatic way once the 
initial conditions are specified. Let us as an example consider the five-line integral MI(5.2) 
(cf. Fig. Hj) which fulfills the following differential equation 

iLM,(5.2) = ^MI(6.2)+ MI(2.2) 

^"^"^^ ■MI(3.4) MI(3.5) 



4m^(z - 1)2 ^ ' ' 8mf{z - 1) 

16mf{z - 4:){z - l)z ^ ' 

(d~ 4:)(z -9)(z + 2) 
^ OA/^ \\Z ,,^ MI(3.10.1) 



16m^{d - 3){z - A){z ~ l)z 

^^'-^^^^ MI(3.10.2)+ ,, MI(4.4). (9) 



Amt{z - A){z - l)'^z 2m4{z 

MI(3.10.1) and MI(3.10.2) denote the Mis of the type 3.10 with a squared massless and 
top quark propagator, respectively. In order to solve this equation it is necessary to know 
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Figure 3: Four- line Mis. For the notation we refer to Fig. [2j 




5.3 5.4 

Figure 4: Five-line Mis. For the notation we refer to Fig. [2l 

the results of all integrals with less than five lines. With the help of the ansatz 

MI(5.2) = 5^Q,fceV(ln^)^ (10) 

the differential equation can be expanded in e and z. As a result it reduces to algebraic 
equations for the coefficients Cj^fc- In every order in e there is one constant Cjjfc which can 
not be determined with this procedure. It it obtained from the initial condition at z = 0, 
which in the case of MI(5.2) can be found in Ref. [35]. In this way we have computed 
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expansion terms up to order z^^ which can be found on the Mathematica file mentioned 
above. For illustration we present the first two expansion terms of MI(5.2) which read 

1 1 In'^ 3 
MI(5.2)mj^ = -7r2ln2-21n^2- -7r2ln3 + 21n221n3-ln21n23 + 



- m — + - In^ 2 
12 2 



,1 1 , 3 , 3 ■ 

in [ In 2 H — In 3 In 2; 

2 4 4 8 



3 , , 3 In^ 3 1 , , In^ 2 , 
+ -In31n2 + - In 2 In 3 In 2 

8 8 4 4 



W 1, 
H L12 , 

48 8 V 4 



0{z'). (11) 



Note that for some integrals the differential equation can be solved with the help of 
Harmonic Polylogarithms [41] which immediately leads to a closed result. 

It is interesting to mention that for the integrals MI(4.3) and MI(4.4) no initial condition 
is needed in order to obtain all the coefficients in the ansatz. They are completely fixed by 
the corresponding differential equation and the solutions for the integrals of the subtopolo- 
gies. For all other integrals initial conditions at z = are required. As already mentioned 
above most of them can be found in the literature or are quite simple to compute using 
standard techniques. However, we could not get analytic results for fiv^ coefficients in 
the e-expansion of the integrals MI(4.5), MI(4.8), MI(5.3) and MI(5.4) at z = 0. We 
calculated these coefficients using the Mellin-Barnes method (see, e.g., Ref [42]) where we 
used the program packages AMBRE [43] and MB [44] . 

The Mellin-Barnes representation for a given integral is not unique. In particular it might 
happen that the convergence of the resulting numerical integration turns out to be good 
in one case whereas a poor convergence is observed in other cases. The crucial quantity in 
this respect is the asymptotic behaviour of the F function for large imaginary part which 
is given by 

6— ►oo 



F(a ± tb) ^ v^e±^5(2-i)e±^^(i°^-i)e-¥6'^-i , (12) 

where the first two exponential factors lead to oscillations. Let us discuss this in more 
detail for the Mellin-Barnes representation of the integral MI(4.5) 

M,(4.5) = ffZVy-p^fe p^,..,....,4-.-» 

\ J J -ZOO 27r2 J_i^ 2m 



2 



^ F(l - e)F(-e - + l)F(-^i) [F(-26 - + 1)] 
^ r(-3e -zi + 2)F(-2e - zi + 2)F(-4e - 2z2 + 2) 
xF(-4e -Z1-Z2 + 2)T{zi - Z2)T{e + Z2)T{2e + Z2) , (13) 



^One more coefficient can be obtained analytically from the requirement that our final result is finite. 
It agrees perfectly with our numerical result. 
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where '~)e = 0.577216 . . . and the contour of integration is chosen in such a way that the 
poles of the F functions with +Zi are separated from the poles of the F functions with —Zi. 
Using the package MB we can expand the integrand in e. For the finite contribution this 
leads to a sum of an analytic part, a one-dimensional Mellin-Barnes integral and a two- 
dimensional one. The latter correspond to the integral in Eq. (|T3l) for e = 0. If we insert 
in this expression the asymptotic behaviour for the F functions as given in Eq. f|T2l) one 
can see that the integrand of the two-dimensional integral falls off exponentially, except 
for Im(z2) = 0,lm(zi) < 0. On this line the drop-off only shows a power-law behaviour 
which is dictated by the last factor of Eq. ( |T2|) . In our particular case the drop-off turns 
out to be extremely slow for the integration contour chosen by MB which corresponds to 
Re(2;i) = —1/4 and Re(2;2) = —1/2. Thus it is hard to get an accurate result by the 
numerical integration since a highly oscillating functions has to be integrated. A closer 
look to the fall-off behaviour in Eq. (fT2!) shows that it is possible to improve the drop- 
off for Im(zi) — > — oo by taking residues of the integrand in Z2 and thus shifting the 
integration contour for Z2 more and more to positive values for Re(z2)- In this way the 
integrand becomes well-behaved and can be integrated numerically with sufficiently high 
precision. 

Let us mention that in the case of MI(4.5) there is an alternative possibility to improve 
the numerical properties of the Eq. (|T3|) : after the variable transformation Z2 ^ Z2 — 
2e MB chooses integration contours which lead to a rapid convergence of the numerical 
integration. We have checked that both approached lead to the same results and obtained 
9 digits for the finite part of MI(4.5) in eleven minutes of CPU time. 

The remaining three integrals show similar properties as MI(4.5). In all cases it is possible 
to end up with integrals which could be integrated numerically. Our results read 



MI(4.5) 



MI(4.8) 



z=0 



m2MI(5.3)|,=o 
m2MI(5.4)| 



2 = 



mt 



2e r 



1 1 
2?^7 



+ 4i7r(l - In 2) 

2 \ 2e 



mt 



1 
2? 



21n2 + i7r -4.81543683(7) 



+ -In21n3 
2 



4 M 4 





19 


237r2 


+ 






y 


24 






- 5 In 2; 




11 


y' 





3 In 3 



2 2 
e 16.690539(1) 



ITtt 



45 3577^ 211n2 2 , , , 

— + + 31n2 + 61n3 + 2Li2(-2) 

2 24 4 ^ ' 



-4Li2 



2.704628(4) - 5.167709(4)i^ 
2.70543(6) - 1.91431(6)i. 



(14) 



The accuracy for the finite part of these integrals is sufficient to obtain the final result 
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with four significant digits. 

Note that contrary to the default settings of MB we do not use Vegas for the multidimen- 
sional numerical integrations. Instead we use Divonne which is available from the Cuba 
library [45]. For the integrals we have considered it leads to more accurate results using 
less CPU time. 

We have performed an independent check of the initial conditions for all the Mis employing 
the method of sector decomposition. In particular we used program FIESTA [46]. 

4 The jtt vertex 

In this Section we discuss 0{aas) corrections to the 7^ vertex due to W boson and gluon 
exchanges with incoming photon momentum at q"^ = 4:mf, the production threshold of 
top quark pairs. This leads to corrections to hf'y mediated by a virtual photon, i.e., to 
the first term of h^/y in Eq.([5]). We denote by the contribution of the sum of all one- 
particle-irreducible diagrams to the 7^ vertex and parameterize the radiative corrections 
in the form 

n = Qt + T'/'''^ + T'A^'''^+T'/'''\ (15) 

where the hat denotes renormalized quantities. Substituting F^ for the Qt in first line of 
Eq. ([5]) and retaining the relevant orders in the electroweak and strong couplings leads to 
the corrections to the helicity amplitudes, hfy\hfy^ and h^^y^ ■ We further decompose 
(and similarly the quantities on the right-hand side of Eq. flTSj) ) according to the 
contributions from the Higgs, Z and W boson exchanges: 

f A = f + f Iz + f • (16) 

In this paper we compute the corrections up to order aag to t^w- The other two 
quantities have been computed in Ref. [8] where the matching coefficients c^'™^ and 
^z,mix j^g^yg Iqqqyl introduced. We have the following relations 

pt,(l,l) _ _««£_^ ^ H,mix f.t,ii,i) _ aos^ Zj^i^ 

^ A,H — 22 ^t'^F^v 5 ^ A,Z — 22 Vt<-^F • U^j 

In our calculation we adapt in the electroweak sector the t'Hooft-Feynman gauge, i.e. 

= I5 which guarantees a simple form of the W boson propagator. Note that 7^ 1 
introduces an additional mass scale in our calculation which would lead to significantly 
more complicated integrals at two-loop order. We want to mention that our final result 
depends on^w This dependence only cancels out after including the self-energy and box 
contributions. There is no gauge parameter dependence in F^ jj which only occurs in the 
vertex and self-energy contributions. Also the vertex corrections involving the Z boson 
(F^ z) are independent of the corresponding gauge parameter. 

Although the one-loop results are well-known, we start our discussion from this order 
since they enter the renormalization of the two-loop expressions. 
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4.1 One- loop corrections 

The renormalized QCD contribution is given by 



,t,(0,l) 



(0,1) 



+ QtZ. 



(0,1) 



where Z2 is the onshell wave function renormahzation for external top quarks. The 
expressions on the right-hand side of Eq. (fT8|) are given by 



r 



t, (0,1) 



^(0,1) 



^* ~A 2 

47r V mf 



CtgCp f jjL 



4:71 \ ml 



t 

2 \ f 



- -4 + e 8 + — 



4-e 

e 



(19) 



where Ci? = 4/3. In Eq. (fT9|) the 0{e) terms are kept since they enter the finite part of 
the two-loop expression. For the renormalized vertex we have the relation 



,t,(0,l) 



Qtc\; 



(1) 



-8g*^(l + 6L,) + 0(e^ 
47r 



(20) 



where = In The one- loop formula for the electroweak corrections is given by 



^t, (1,0) 
A 



pt, (1,0) ^ ryK^^V) rji-6 



(1,0) 



n3 v(l,0) 



(21) 



where .^^^^ ^ counterterm associated with Z-photon mixing at zero-momentum trans- 
fer [47] which reads 



7(i>o) 



a 



1 , 



w 



+ 0{e) . 



(22) 



As we will see later the 0{e) term for Z^rp^ is not needed. We present the remaining 
two ingredients as a series expansion of z = M^/mf, the exact formulae are collected in 
Appendix for convenience. The results for Z^^^ reads 



7(1,0) 



a f jj 



+ e 



,2 \ ^ 



96z 8 











in 




A 






si " 


4 


+ 


57r2 
16 


T 


3 

^ 4 



iTT 1 

y ~ 4 



in 21 

y ^ 16 



+ 



3L, 



TT 

+ T7T - 



QzTT 27 



16 16 12 8 



32 



(23) 
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and is given by 



■^t, (1,0) 
AW 



a 



Qb ( 32 



mt 



1/1 1\ Ifin 
7 81 + 4 +1'- 



21n2 1 

3 ^ 6 



in 1 

^ T ^ Y ~ 12 



21n2 



21n2 



7 

48 



+ e 



l/77r2 i7r(41n2-l) 21n2(31n2-5) 5 



z 



96 



+ 



6 



9 



— + — 

8 4 

16 16 



16 



m(41n2-l) 21n2(31n2 + l) 7 
6 ^ 9 ^ 36 



21n2(31n2 + l) vr^ 
9 ^ 12 



(I61n2-l) 115 



+ 



a 



rp 3 

47rs2 * 



1 



+ 



2; 



T 



24 

21n2 2 

— + 3 



288 



+ 



31n2 + 3 



5L, 



15 5i7r 

16 ^ ~8~ 



71n2 



+ e 



1 /77r2 2i7r(ln2-l) 21n2(31n2-8) 13 
i lis" ^ 3 9 Y 



77r2 



-3i7r(ln2- 1) +31n2(ln2-2) +6 



+ 



5L? 



IIL^ 57r2 



16 



16 



+ 



12 



+ 



17T 



(71n2 + l) 71n2 2 39\ 



H ] z 

32 J 



(24) 



where = Inz and a = e^/(47r) is the fine-structure constant in the Thomson hmit. 
Again the 0{e) terms are retained due to their relevance for the two- loop renormalization. 

^ ,y diiLi ■ -^^ one-loop order this 



Note that we keep the imaginary parts of both F^' w^'* and zi^'^^ 



only affects the finite part; at two loops also the pole parts are concerned (see below). 

We refrain from listing an analytical result for f ^' ^^'^^ but compare in Fig. \5\ the approx- 
imated result to the exact one. The latter is represented by the (red) solid line whereas 
the dashed lines correspond to the expansions including successively higher orders in z. 
As one can see, the expression including the correction of order z^ provides at the physical 
point z ~ 0.23 a perfect approximation to the exact result far below the per cent level. 
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Figure 5: Electroweak one-loop corrections to the photon vertex F^^^"-* normahzed to 
a/ {Anslj) / z as a function oi z = M^/rn^, for real part (left panel) and imaginary part 
(right panel), respectively. The (black) dashed lines include successively higher orders in 
z starting from (long dashes) to z^ (short dashes) and the (red) solid line is the exact 
one-loop correction. The (orange) band marks the physical range oi z = (80.40 GeV/rrit)'^ 
varying the top quark mass between 165 and 175 GeV. 



The approximated results are based on the following expressions 



.t,{i,o) 

AW 



a 



AttsI 



0.20 



+ (0.48 + 0.79 z + 0.25 In^) 



+ z{- 0.0024 - 1.37 i - 0.44 lnz)+z'^{- 0.072 + 1.39 i + 0.44 In z) 
+ z^ (0.34 - 0.53^-0.17 In + z'^ (0.13 + 0.23 i + 0.074 In^) 

+ 0{z^) 

2.48i/^ + 0.25i + 0.39^ - 0.095^2 + 0.017^3 + 0.00013^4 



+ i (2. Ill - 0.80^ + 0.18^2 - 0.015^3 + 0.0014^*) 



X 10" 



(25) 



where the subscript in the last line indicate their order in the z expansion and for the 
input parameters the following values have been used [48,49] 



a, = 0.108, a(Mz) = 1/128.9, 
Mw = 80.40 GeV , rrit = 172.4 GeV . 



0.23. 



(26) 



Note that in our numerical analysis we use a at high energy scalejfl 

The corrections in Eq. (!25|) are dominated by the leading terms proportional to m^/M^. 
One observes a rapid convergence, so that the term of order z^ can safely be neglected. 
Inserting the results in Eq. ( JT5l) the overall size of the electroweak corrections (from the 
diagrams involving a W boson) amounts to about 0.5% and is unusually small. For 



^This is theoretically preferable because it is devoid of non-perturbative hadronic effects. 
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comparison, we note that f and f^'^]^°^ lead to corrections of 0.3% and 3.2% (for 
Mh = 120 GeV), respectively. Let us mention that the one-loop QCD corrections provides 
a contribution "—61 x 10~^" to the last line of Eq. fl25|) thus resulting in a 9% correction. 

From Eq. fl25p one obtains the corresponding corrections to the helicity amplitude as 

= afl'fi)°\ (27) 

which immediately leads to the correction to the cross section with the help of Eq. ([T]). 
Taking at tree-level both the photon and Z exchange diagram we obtain a shift of 0.9% 
to -R(e^e~ tiX) due to the W boson contribution to 7tt vertex at one-loop. 



4.2 Two-loop order aas renormalization 

The renormalized 7tt vertex at order aas is given by 

_ 3 (1,1) ^ (1,1) 

+ r*'(°-^) + V'f^'A + % f*'(°'^\ (28) 

where the first line corresponds to genuine two-loop diagrams and the second line consists 
of products of one-loop diagrams. For the latter we already listed all the relevant expres- 
sions in the previous Subsection. Note that in the last term the renormalized one-loop 
vertex appears and thus the 0[e) term for Z^^^^ is not needed. 

Formula fl28|l takes only care of the renormalization of the external lines and the elec- 
tric charge which means that the un-renormalized two-loop quantities in the first line 
are understood as the sum of the amputated two-loop diagrams and the corresponding 
counterterm diagrams for the top quark mass and the top quark Yukawa coupling. The 
latter are renormalized in the onshell scheme. 

It is easy to see that for the two-loop counterterm we have Z^^'^^ = since at one-loop order 
only bosonic and no fermionic diagrams contribute. The two-loop onshell wave function 
factor Z2^'^^ has been computed in Ref. [50]. We confirmed the result by an independent 
calculation and added the imaginary part which is necessary in our framework. The result 
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reads 



7(1.1) 
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(29) 



where terms up to order z^ have been included (C(3) = 1.20205 ■ ■ ■ )• 

In the following we provide the result for the un-renormalized vertex corrections where 
the finite part is given in numerical form. Our result reads 



A,W 



a agCp [2 1 



47rs^ 47r 
17L 
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16 
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96 



48 



+ 2-66 - 3.79] - [1.13 - 1.57i] ) + 4.00 

+ - 0.50L2 - [1.98 - ll.OOi]) + 0.13L2 + 0.08L^ + [10.03 - 23.63z] 

(^[5.24- 6.68i] -2.13L^) +0.53L^-L^ [0.77- 7.76i] 



+ z 



[11.84 + 5.62i] 



+ z' 



(^[1.23 - l.lli] -0.35L, 



+ 0.09 Ll - L, [3.80 + 1.16i] + [9.31 + 18.63i] 



+ 0{z' 



(30) 



Let us mention that in our calculation we allowed for a general QCD gauge parameter 
^5 and used the independence of r^'*'^'^'' as a welcome check for the correctness of our 
result. Note that for the cancellation of ^5 it is important to include the counterterm 
diagram for the top quark mass. The remaining ingredients in Eq. (128|) are individually 
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^s'-independent. A further check of our calculation is based on a setup where we choose 
Mw = from the very beginning. This leads to significantly simpler expressions during 
the reduction to master integral, which is completely independent from the one for finite 
Mw 



4.3 0{aas) corrections to the ^tt vertex 

Inserting all ingredients into Eq. (l28l) leads to 

ft,{l,l) _ " CXsCf 



AW 



47rs^ 4:71 



0.45 -i 2.06 



6.34-^25.14- 2.00 Inz 



+ z 1^-6.27-^6.10+ (2.16 + i4.10) In^^ 

+ 1^13.50 - i 31.29 - (4.53 + i 2.9l) In 

+ (^-41.06-^5.55 + (0.59 + i 13.10) \nz^ 

+ / (17.11 -i9.99 - (l.86 + i5.36) \nz] + 0{z^ 



( - 0.64i/^ + 2.89i - 0.64^ + 0.30^2 - 0.13^3 + 0.014^4) 
+ 2.91i/^ - 7.73i - 0.83^ - 0.39^2 - 0.081^3 - 0.0013^4) 



X 10" 



(31) 



In Fig.[6]we show the result for f^*^^'^^ including first five terms of the z expansion. Taking 
the difference of two successive curves as a measure for the quality of the approximation 
we observe a rapid convergence at the physical point. Note that in contrast to the one- loop 
case the leading 1/z contribution is numerically not dominant. 

In analogy to Eq. (|27l) we obtain for the correction to the helicity amplitude 



A,W 



'''ly 



which results in a correction of 0.1% to Rln 



Ve-L A,W 



(32) 



Let us in the following briefly compare the new vertex corrections to the ones induced 
by a Higgs and Z boson. Note that the latter two contain a non-trivial scale dependence 
which is canceled by the corresponding contribution from the effective theory [8]. Choosing 



^We did not take into account 0{aas) interference terms, e.g., 5R ^ (Stt/s) [2/i^^y'ft,j''y''iJy] . Such 
correction should be considered once the the two-loop box contributions are available. 
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0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

z z 



Figure 6: Electroweak two-loop correction to the photon vertex T^/^l(^^ normahzed to 
{agaC f) / {IGii^ sIj) / z as a function of z = M^/m^, for real part (left panel) and imaginary 
part (right panel), respectively. The (black) long-dashed lines include successively higher 
orders in z starting from z° (long dashes) to z^ (short dashes) and the orange band is the 
physical range of z (see Fig. E]). 

fj, = rrit one obtains [6,8] 

21.1 X 10"^ (10.6 X 10~^) for Mh = 120 (200) GeV, 
3.0 X 10^^ 
1.7 X 10-^ 

-17.6 X 10~^ (-6.6 X 10"^) for Mh = 120 (200) GeV, 
0.2 X 10"^ 

-1.0 X 10"\ (33) 

One observes quite small corrections from the W and Z boson induced contributions. 
From Eq. (|33il one can read off that relatively big one-loop effects are obtained for light 
Higgs boson masses. However, there is a strong cancellation between the one- and two- 
loop terms resulting in corrections which have the same size as the sum of the one- and 
two-loop contributions of the W and Z boson diagrams. In general moderate effects are 
observed suggesting that in the electroweak sector perturbation theory works well, which 
is in contrast to the pure QCD corrections. 

Let us at this point comment on the imaginary parts contained in Eqs. ( l25l) and ( l3Ti) 
which are not taken into account in the numerical estimates for the corrections to R 
presented above. As we mentioned previously, "Im" in Eq. applies to the imaginary 
part which corresponds to the tiX final state or experimentally indistinguishable cuts 
involving bottom quarks and W bosons. Thus, it is necessary to separate imaginary parts 
arising from cutting, e.g., two W boson or two b quark lines from the tiX cuts in order 
to make a phenomenological prediction. This requires a dedicated analysis in the loop 
calculation, which is beyond the scope of this paper (see, e.g., Ref. [11]). 



A,H 

t,(lfi) 
A,W 

A,Z 

A,W 
A,Z 
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5 Conclusions and outlook 



Mixed two-loop electroweak/QCD corrections to the 7tt vertex due to W boson and 
gluon exchange have been computed. The new contribution completes the order aOg 
corrections to the 'jtt vertex. The numerical evaluation leads to a shift of 0.1% in the 
threshold production cross section of top quark pairs at e'^e~ colliders which is small 
as compared to the aimed 3% uncertainty for the theory predictions. Nevertheless, it is 
remarkable that in the sum of the order a and order aa^ correction terms the sizeable 
one-loop contribution from the Higgs boson induced diagrams is screened resulting in 
numerical values comparable to the W and Z boson contributions. 

We want to mention that the corrections evaluated in this paper can be taken over in a 
straightforward way to the vector coupling of the Zti vertex. Note, that the axial- vector 
contribution is suppressed at threshold and thus one- loop corrections are sufficient. 

The only missing building block in order to complete the order aag corrections to the 
process e'^e~ — > tt are the two-loop box diagrams. They are technically more involved 
and are thus postponed to future work. 
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Appendix: Exact result for 



t,(i,o) 



w 



Keeping the full dependence on e the exact one-loop result for the ■ytt vertex due to the 
W boson exchange reads (for massless bottom quarks) 



-.4,(1,0) 
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Ansl V -3 + 2 e 
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(l-^)(l + 2(l-e)^) „ „ ^ 



8z{l + z) 



«,(i±l(izi)|i^B<«)(4™?.0.0) 



(34) 
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where Qb = —1/3 is bottom quark electric charge normahzed to the one of positron. The 
corresponding contribution to the wave function renormahzation constant reads 



,(1,0) 



a 



47rs? 



(1 + 2(1-6) z) 



A^o\m^) (l + ;,)(l + 2(l-e)z) 
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(35) 



The loop-functions Aq'' and Bq^^ are given by 
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(36) 
(37) 



(3^ 



(39) 

where an analytic continuation by (m^ + iO) is understood. Expansions with respect to 
e of the Gauss-hypergeometric functions 2-^1 around integer values is well known in the 
literature (see, e.g., Ref. [38] or the package HypExp [51]). 
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